Electron-phonon interaction in the lamellar cobaltate Na^CoC^ 
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We study theoretically and experimentally the dependence of the electron-phonon interaction in 
Na^CoC^ on the sodium concentration x. For the two oxygen phonon modes found in Raman 
experiments, A\ g and Ei g , we calculate the matrix elements of the electron-phonon interaction. 
Analyzing the feedback effect of the conduction electrons on the phonon frequency we compare the 
calculated and experimentally observed doping dependence of the Ai g mode. Furthermore, due 
to the momentum dependence of the electron-phonon coupling for the E\ g symmetry we find no 
renormalization of the corresponding phonon frequency which agrees with experiment. Our results 
shed light on the possible importance of the electron-phonon interaction in this system. 
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Introduction. The origin of the unconventional su- 
perconductivity in low-dimensional perovskite systems 
attracts much attention and belongs to the most chal- 
lenging questions of condensed matter physics. The 
best known example among these materials are high-T c 
cupratc superconductors. There, one of the possible sce- 
narios for the Cooper-pair formation is the so-called spin 
fluctuation mechanism. At the same time, due to the 
complexity of the transition metal oxides, other energy 
scales arc present, and their role in the formation of su- 
perconductivity remains under debate. This, in particu- 
lar, concerns the electron-phonon interaction. For exam- 
ple its relevance for superconductivity in layered cuprates 
and the anomalous normal state has been discussed in, 
e.g., Ref. 0. Despite some progress, a complete under- 
standing of the physics of electron-phonon coupling in 
perovskites is still lacking because of the crystallographic 
complexities of these materials. 

The discovery of superconductivity with T C =4.6K in 
water intercalated sodium cobaltate, Na x Co02 • 2/H2O 
Q, is of great interest on its own and also because of 
similarities with layered cuprates. The sodium cobaltate 
has a quasi-two-dimensional layered structure with C0O2 
layers and rich phase diagram as a function of the Na con- 
centration, which includes superconductivity at x « 0.3, 
an insulating phase at x ~ 0.5, and unusual magnetism 
for x > 0.6 [3j. There is also increasing experimental and 
theoretical evidence for unexpected strong correlation ef- 
fects as the cobaltates approach the band insulting limit 
at x = 1 @,0,0,0,0,0,l!flElL In Na ^ Co °2 the Na ions 
reside between the C0O2 layers, with Co ions forming a 
triangular lattice, and donate x electrons to the partially 



filled Co-d(t2 g ) orbitals. Due to the presence of a trigonal 
crystalline electric field, the t2 g levels split into the higher 
lying a\ g singlet and the two lower lying e' g states [12| |. 
Angle- Resolved Photo-Emission Spectroscopy (ARPES) 
reveals a doping dependent evolution of the Fermi 



surface, which shows no sign of the e' g hole pockets for 
0.3 < x < 0.8. The observed Fermi surface is centered 
around the T point and has mostly a\ g character. It has 
been argued that such an effect may arise due to strong 
electronic correlations [E EH , or Na induced disorder , 
however, no consensus in the literature has been reached 
yet (sec, for example, [IE EE EH ) • 

Despite of intensive studies of the electronic and mag- 
netic properties little is known about the phonon ex- 
citations and their doping evolution in Naa;Co02. At 
the same time, due to the relatively low superconducting 
transition temperature the possible relevance of phonons 
for superconductivity cannot be neglected. For exam- 
ple, the role of the electron-phonon coupling in Na K Co02 
has been discussed in the context of its relevance to su- 
perconductivity and charge ordering on the basis of a 
t — V model [201 ] . In addition, due to some similarity 
with high-T c cuprates the understanding of the phonon 
renormalization in the sodium cobaltates is of great im- 
portance. Initially, the effect of renormalization of the 
optical phonons by the conduction electrons in layered 



cuprates has been considered in Ref. [21 



In this Rapid Communication we investigate the 
electron-phonon interaction in the Na K Co02 as a func- 
tion of doping concentration and its superconducting rel- 
ative by means of Raman spectroscopy. We observe two 
oxygen phonon modes at small wave vectors with A\ g 
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and E\ g symmetries. Then we derive the diagonal and 
off-diagonal electron-phonon matrix elements for these 
modes. Calculating the renormalization of the phonon 
frequencies by conduction electrons we compare our re- 
sults with the doping dependent evolution of the A\ g 
mode and obtain the electron-phonon coupling constant 
g a j 9 f = 3mcV. Due to the structure of the electron- 
phonon matrix element for the E\ g mode we obtain 
no doping dependence of the corresponding phonon fre- 
quency in good agreement with experiment. Our results 
shed light on the possible role of the electron-phonon in- 
teraction in this compound. 

Experiment details. Raman scattering experiments 
have been performed in quasi-backscattering geometry on 
freshly cleaved single crystal surfaces. The sample have 
been fully characterized using basic thermodynamic as 
well as spectroscopic techniques [22], [H, [24, 25, 2|| 27 1. 
After cleavage the crystals were rapidly cooled down in 
He exchange gas to prevent degradation. In Na^CoC^ in- 
plane E\ g and out-of-plane A\ g ox yge n modes have been 
observed in Raman scattering 

HE 01 and the corre- 
sponding oxygen displacements are depicted in Fig. [IJa) 
and (b). The Co site is not Raman-active. Modes of 
the Na sites have not been identified unambiguously . 
This is probably related to disorder on the partially filled 
sites. The two-dimensionality with respect to structure 
and bonding leads to a decoupling of the Na and the 
C0O2 layers. The observed doping-dcpendcncc of the 
A\ g and E\ g oxygen phonon frequencies are shown in 
Fig. [2j The cross-over from one to the other crystallo- 
graphic phases (shaded areas) given by a different occu- 
pation of the Na sites leads for the A lg modes to small 
additional frequency shifts and for the E lg modes to new 
modes with a larger energy off-set. The latter are omitted 
for clarity. The two phonon modes display a markedly 
different doping dependence. 

Tight-binding model. To describe the electronic sub- 
system we use a tight-binding <2g-band model with pa- 
rameters (in-plane hoppings and the single-electron en- 
ergies) derived previously from the ab-initio LDA (Lo- 
cal Density Approximation) calculations using projection 
procedure for x = 0.33 Q. 

The free-electron Hamiltonian of the i2g-band model 
in a hole representation is given by 



H 



E 



k.cr Q,/3 



where n kQ(T = d^ aa dv. aa , d^ aa (dl aa ) is the annihilation 
(creation) operator for the t2g-hole with spin er, orbital 
index a, and momentum k, t^f is the hopping matrix 
element, and e a is the single-electron energy. To obtain 
the dispersion we diagonalize the Hamiltonian ([I]) calcu- 
lating the chemical potential /i self-consistently. Due to 
the non-zero inter-orbital hopping matrix elements. a\ g 
and e' g bands are hybridized. However, only one of the 
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FIG. 1: (Color online) Schematic illustration of the A± g (a) 
and Eig (b) phonon modes. The arrows indicate the oxygen 
displacement in the CoOe octahedra. On the left side of (a) 
we indicate the crystallographic a, b, and c directions, (c)-(d) 
The calculated momentum dependence of the electron-phonon 

structure factors, [fq ] , in the first BZ for the A\ a and E\ g 
phonon modes, respectively. 



hybridized bands crosses the Fermi level. We refer to the 
diagonalizcd bands as with the new orbital index a' . 

Electron-phonon interaction. In analogy to previ- 
ous considerations for cuprates (3C" 



31|, we derive the 



electron-phonon matrix elements, <7kq, for the Ai g and 
Ei g phonon modes depicted in Fig.QJa) and (b). Namely, 
to obtain the main contribution to the diagonal (intra- 
band) part of the electron-phonon interaction we ex- 
pand the Coulomb energy between Co and oxygen 
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the small displacements of the oxygen ions. Here, e is the 
electron charge, e* = — 2e is the oxygen ion charge, e is 
the dielectric constant, Ri are the Co ion positions, rj j7 
are the vector positions of the vibrating oxygens, and in- 
dex 7 (—7) labels the three oxygen positions within CoOg 



unit cell above (below) the Co layer. Here, 



refers 



to the diagonal form of the Hamiltonian fl}. After intro- 
ducing the creation (annihilation) operator 6_ q (6 q ) for 
the phonon with momentum q, we arrive to the following 
form of the electron-phonon interaction 



H 



diag 



ei- ph = E ^ V cLvCk- q «v(& q + &L q ). (2) 

k,q,a',<r 

For the sake of simplicity we assume the diagonal 
electron-phonon interaction is independent on the orbital 
index a'. Thus, for the Ai g and E\ g optical Raman- 
active phonon oxygen modes one finds c/ q 1 
.9q ls = g^-Fa 19 1 where the structure factors of the 
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electron-phonon interaction are 



F> = cos 91 " 92 ' — 91 + 93 



600 



Fa 



3 

gl - 92 _ 1 

3 2 



3 

qi + ?3 



92 + 93 

cos , (3) 



3 ' 

92 + 93 



(J) 



Here, qi = (\/%/2)q x -\q v , q 2 = q y , 93 = (\/%/2)q x + \q y , 
in units of 27r/a with a being the in-plane lattice con- 

I— — where lo t is the 



Stant, = ~^£TfjWy 2Mu r 

corresponding bare phonon frequency (r = ^4i ff ,£i s ), 
L/A lg = d = a/v6 is the distance between the Co and the 
oxygen plane, Le 19 = I = a/V3 is the planar distance be- 
tween Co and oxygen, and M is the oxygen mass. Assum- 
ing that in the band insulator, Na a;= iCo02, the renor- 
malization of the phonons by the conduction electrons is 
absent, we use loa 1 = 589cm _1 and u>e 1 — 470cm _1 . 
These values are close to those obtained by the first prin- 
ciples calculations [32[ ■ The resulting momentum depen- 
dence of the structure factors for the both modes is shown 
in Fig. QJc) and (d). Interestingly, one sees that while 
the <7 q for the Ai g mode shows a maximum at the BZ 
center, the corresponding g q for the E\ g mode vanishes 
there. Therefore, for q = the electron-phonon coupling 
for the E\ g channel is zero. Taking e ~ 20 we estimate 

9t: 9 * 0.05eV. 

The off-diagonal (interband) contribution to the 
electron-phonon interaction arises mainly from the mod- 
ulation of the inter-orbital Co-Co hopping matrix ele- 
ment via oxygen. Assuming the linear terms in the ex- 
pansion of the nearest neighbors hopping matrix element 
t"j(u 7 ) = iy + V a/3 u 7 over the oxygen displacements 
Uj, one obtains: 



H 



off \p 

el—ph / j 
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L'<r C k-q/3' CT (6q + &l q ). (5) 



where g^' = ( 7 (k) + 7 (k + q)) with 7 (k) = 

cos &2 + cos ^3 + cos fci being the Co lattice structure fac- 
tor. Again one could see that for q = the off-diagonal 
electron-phonon coupling for the E\ g channel is zero due 
to the momentum dependence of the electron-phonon 
structure factor, Fq l9 . Therefore, in Raman experiments 
which probes q — response this mode shows no doping 
dependence due to the coupling to the electronic subsys- 
tem. This is also confirmed by the fact that the observed 
phonon mode energy for all doping levels lies close to the 
value obtained by ab-initio calculations 32J. The only 
Raman-active optical phonon mode which will couple to 
the conduction electrons at q = is the A\ g mode. 

In the following we consider the renormalization of the 
A\ g phonon. The corresponding Dyson equation reads: 



D- 1 (q,(j) = Do 1 (Lj)-U(q,uj) 



(0) 



where Dq(u>) 



2oj r 



u 2 -Lu r2 +i5 * s ^ ne m omentum-independent 
bare phonon propagator. The polarization operator is 
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FIG. 2: (Color online) The experimental Raman data at 10K 
is shown by (blue) circles. The calculated doping dependence 
of the A\ g (a) and E\ g (b) phonon modes is shown by red 
crosses, the solid curve is a guide to the eye. The slight scat- 
tering of the experimental points around x — 0.5, x — 0.76 
and x = 0.9 (shown by the dashed areas) are the result of 
the different crystallographic phases around these points [27j|. 
These structural modifications only weakly affect the electron- 
phonon coupling as it mainly involves Na ordering. Measure- 
ments of Ai g mode at 290K (not shown) follow a similar trend 
as the 10K data. 



given by: 



n(q,w) = -2 



EE 

a',/3' k 



a a' 



-q)-/0£) 



lu — e 



k+q 



(7) 

where /(e) is the Fermi function. To find the renormal- 
ization of the bare phonon frequency and to compare the 
results to the Raman experiments we set q — > limit 
and solve Eq. ^ numerically as a function of the dop- 
ing concentration. The main contribution to the renor- 
malization of the optical phonon modes comes from the 
interband transitions, i.e. terms with <7k q while intra- 
band transitions renormalize the acoustic phonon modes. 
The results of our numerical calculations are shown in 
Fig. [21 The doping evolution has been deduced by cal- 
culating II(q — > 0, a;) for various x values. We obtain 
the value of the off-diagonal electron-phonon interaction 
g ff ~ 3meV by comparison to the experimental data 
points. This value is an order of magnitude smaller than 
the diagonal contribution to the electron-phonon interac- 
tion. One sees that the phonon renormalization changes 
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smoothly as a function of doping. The rcnormalization 
effects tend to vanish close to the band insulator regime 
at x = 1 because all of the Fermi functions entering 
Eq.© are equal one and ReII(q, 0) = 0. Away from 
this point Rell will be simply determined by the number 
of holes in the system. This explains the linear doping 
dependence of the renormalization effects. We also note 
that the electronic correlation effects seem to play a mi- 
nor role for the renormalization of the Raman-active op- 
tical phonon modes. The situation may differ, however, 
for the acoustic phonons where the density of the states 
at the Fermi level plays the most important role (33j . 

Of course, our estimated value for the electron-phonon 
interaction corresponds to the intcrband scattering which 
is not directly related to the superconductivity. There- 
fore, we use our estimated value for the diagonal part 
of the electron-phonon interaction and obtain using the 
BCS formula k B T c = l.Uoj Alg exp {-\/N(EF)g M ») 
T c ~ 7 K for the N{E F ) w 4.0" states/eV 0. This 
estimate is in agreement with the observed T c in water 
intercalated Na^CoC^, which points towards potential 
relevance of the electron-phonon interaction for the su- 
perconductivity in this compound. At the same time, the 
smooth evolution of the phonon frequencies as a function 
of the doping concentration on the one hand, and the rich 
phase diagram of the lamellar cobaltatcs, on the other, 
requires further understanding the role of the phonons 
in the formation of superconductivity One interesting 
possibility to explore is the possible enhancement of T c 
due to resonance levels the may be introduced by the 
water molecules [34, 3^, 



36| . Very recently the so called 
Suhl-Kondo resonance was suggested as an origin of the 
Cooper paring in the Na x Co02 • J/H2O [33] . In any case 
a study of the isotope effect on T c is highly desirable. 

We would like to thank P. Fuldc, B. Keimer, I. Mazin, 
K. Morawetz, E. Schneider, A. Yaresko, and G. Zwick- 
nagl for useful discussions. We thank P. Scheib and A. 
Boothroyd for contributions to our experimental study. 
M.M.K. acknowledges support from INTAS (YS Grant 
05-109-4891) and RFBR (Grants 06-02-16100, 06-02- 
90537-BNTS). I.E. acknowledges support from Volkswa- 
gen Foundation. The experimental studies have been 
supported by DFG and ESF-HFM. 



[1] A. Lanzara, P.V. Bogdanov, X.J. Zhou, S.A. Kellar, D.L. 

Feng, E.D. Lu, T. Yoshida, H. Eisaki, A. Fujimori, K. 

Kishio, J. -I. Shimoyama, T. Noda, S. Uchida, Z. Hussain, 

and Z.-X. Shen, Nature 412, 510 (2001). 
[2] K. Takada, H. Sakurai, E. Takayama-Muromachi, F. 

Izumi, R.A. Dilanian, and T. Sasaki, Nature 422, 53 

(2003). 

[3] M.L. Foo, Y. Wang, S. Watauchi, H.W. Zandbergen, 
T. He, R.J. Cava, and N.P. Ong, Phys. Rev. Lett. 92, 
247001 (2004). 



[4] 

[5] 

[6 
[7 
[§ 

[10 

[11 

[12 
[13 

[14 

[15 
[16 

[17 

[18 

[19 
[20 



[21 
[22 



[23 
[24 

[25 
[26 

[27 

[28 



C. Bernhard, Ch. Niedermayer, A. Drew, G. Khaliullin, 
S. Bayrakci, J. Strempfer, R.K. Kremer, D.P Chen, C.T. 
Lin, and B. Keimer, Europhys. Lett. 80, 27005 (2007). 
M. Lee, L. Viciu, L. Li, Y. Wang, M.L. Foo, S. Watauchi, 
R.A. Pascal, Jr., R.J. Cava, and N.P. Ong, Nat. Mater. 
5, 537 (2006). 

J. Merino, B.J. Powell, and R. H. McKenzie, Phys. Rev. 
B 73, 235107 (2006). 

J.O. Haerter, M.R. Peterson, and B.S. Shastry, Phys. 
Rev. Lett. 97, 226402 (2006). 

M.M. Korshunov, I. Eremin, A. Shorikov, V. I. Anisimov, 
M. Renner, and W. Brenig, Phys. Rev. B 75, 094511 
(2007). 

M. Daghofer, P. Horsch, and G. Khaliullin, Phys. Rev. 
Lett. 96, 216404 (2006). 

C. A. Marianetti and G. Kotliar, Phys. Rev. Lett. 98, 
176405 (2007). 

Meng Gao, Sen Zhou, and Ziqiang Wang, Phys. Rev. B 
76, 180402(R) (2007). 

D. J. Singh, Phys. Rev. B 61, 13397 (2000). 

M.Z. Hasan, Y.-D. Chuang, D. Qian, Y.W. Li, Y. Kong, 
A. P. Kuprin, A.V. Fedorov, R. Kimmerling, E. Roten- 
berg, K. Rossnagel, Z. Hussain, H. Koh, N.S. Rogado, 
M.L. Foo, and R.J. Cava, Phys. Rev. Lett. 92, 246402 
(2004). 

H.-B. Yang, S.-C. Wang, A.K.P. Sekharan, H. Matsui, S. 
Souma, T. Sato, T. Takahashi, T. Takeuchi, J.C. Cam- 
puzano, R. Jin, B.C. Sales, D. Mandrus, Z. Wang, and 
H. Ding, Phys. Rev. Lett. 92, 246403 (2004). 
S. Zhou, M. Gao, H. Ding, P.A. Lee, and Z. Wang, Phys. 
Rev. Lett. 94, 206401 (2005). 

D.J. Singh and D. Kasinathan, Phys. Rev. Lett. 97, 
016404 (2006). 

H. Ishida, M.D. Johannes, and A. Liebsch, Phys. Rev. 
Lett. 94, 196401 (2005). 

L. Balicas, J. G. Analytis, Y. J. Jo, K. Storr, H. Zand- 
bergen, Y. Xin, N. E. Hussey, F. C. Chou, and P. A. Lee, 
Phys. Rev. Lett. 97, 126401 (2006). 

G. J. Shu, Andrea Prodi, S. Y. Chu, Y. S. Lee, H. S. 
Sheu, and F. C. Chou, Phys. Rev. B 76, 184115 (2007). 
A. Foussats, A. Greco, M. Bejas, and A. Muramatsu, 
Phys. Rev. B 72, 020504(R) (2005); M. Bejas, A. Greco, 
and A. Foussats, Phys. Rev. B 75, 033101 (2007). 
R. Zeyher and G. Zwicknagl, Z. Phys. B 78, 175 (1990). 
P. Lemmens, V. Gnezdilov, N.N. Kovaleva, K.Y. Choi, H. 
Sakurai, E. Takayama-Muromachi, K. Takada, T. Sasaki, 
F.C. Chou, D.P. Chen, C.T. Lin, and B. Keimer, J.Phys.: 
Condens. Matter 16, S857 (2004). 

P. Lemmens, K.Y. Choi, V. Gnezdilov, E.Ya. Sherman, 
D.P. Chen, C.T. Lin, F.C. Chou, and B. Keimer, Phys. 
Rev. Lett. 96, 167204 (2006). 

S. P. Bayrakci, C. Bernhard, D. P. Chen, B. Keimer, R.K. 

Kremer, P. Lemmens, C. T. Lin, C. Niedermayer, and J. 

Strempfer, Phys. Rev. B 69, 100410(R) (2004). 

D.P. Chen, H.C. Chen, A. Maljuk, A. Kulakov, H. Zhang, 

P. Lemmens, and C.T. Lin, Phys. Rev. B 70, 024506 

(2004). 

C.T. Lin, D.P. Chen, P. Lemmens, X.N. Chang, A. 
Maljuk, and P.X. Zhang, Journ. of Cryst. Growth 275, 
606 (2005). 

C.T. Lin, D.P. Chen, A. Maljuk, and P. Lemmens, Journ. 
of Cryst. Growth. 292, 422 (2006). 

M.N. Iliev, A. P. Litvinchuk, R.L. Meng, Y.Y. Sun, J. 
Cmaidalka, and C.W. Chu, Physica C 402, 239 (2004). 



5 



[29] P. Lemmens, P. Scheib, Y. Krockenberger, L. Alff, F.C. 

Chou, C.T. Lin, H.-U. Habermeier, and B. Keimer, Phys. 

Rev. B 75, 106501 (2007). 
[30] A. Nazarenko and E. Dagotto, Phys. Rev. B 53, R2987 

(1996). 

[31] S.G. Ovchinnikov, and E.I. Schneider, Zhur. Eksp. Teor. 

Fiz. 128, 974 (2005) [JETP 101, 844 (2005)]. 
[32] Z. Li, J. Yang, J.G. Hou, and Q. Zhu, Phys. Rev. B 70, 

144518 (2004). 

[33] O.V. Dolgov, O.K. Andersen, and I.I. Mazin, 
larXiv:0710.066TV l (unpublished). 



[34] W. A. Little, Phys. Rev. 134, 1416 (1964). 

[35] V. L. Ginzburg, Zh. Eksp. Teor. Fiz. 47, 2318 (1964); 
[Sov. Phys.-JETP 20, 1549 (1965)]. 

[36] V. L. Ginzburg, Contemp. Phys. 9, 355 (1968) [INSPEC]; 
Usp. Fiz. Nauk 101, 185 (1970) [CAS]; [Sov. Phys.-Usp. 
13, 335 (1970)]; Ann. Rev. Mater. Sci. 2, 663 (1972); 
V. L. Ginzburg and D. A. Kirzhnits, Phys. Rep. 4, 345 
(1972). 

[37] K.Yada and H.Kontani, larXiv:0801.3495V l (unpub- 
lished); J. Phys. Soc. Jpn. 75, 033705 (2006). 



